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ABSTRACT 



Context. Helically forced magneto-hydrodynamic shearing-sheet turbulence can support different large-scale dynamo modes, al- 
though the a -Si mode is generally expected to dominate because it is the fastest growing. In an a -CI dynamo, most of the field 
amplification is produced by the shear As differential rotation is an ubiquitous source of shear in astrophysics, such dynamos are 
believed to be the source of most astrophysical large-scale magnetic fields. 

Aims. We study the stability of oscillatory migratory a -SI type dynamos in turbulence simulations. 

Methods. We use shearing-sheet simulations of hydromagnetic turbulence that is helically forced at a wavenumber that is about three 
times larger than the lowest wavenumber in the domain so that both a -CI and dynamo action is possible. 

Results. After initial dominance and saturation, the a -SI mode is found to be destroyed by an orthogonal mode sustained by the 
helical turbulence alone. We show that there are at least two processes through which this transition can occur 
Conclusions. The fratricide of o- - f2 dynamos by its sibling is discussed in the context of grand minima of solar and stellar activity. 
However, the genesis of o- - Q dynamos from an a- dynamo has not yet been found. 

Key words. Sun: dynamo - magnetohydrodynamics (MHD) 



1. Introduction 

The observed existence of large-scale astrophysical magnetic 
fields, for example galactic or solar fields, is usually explained 
by self-excited dynamo action within electrically conducting flu- 
ids or plasmas. However, this mechanism of field amplification 
continues to be a matter of debate as the existing theory encoun- 
ters problems when extrapolated to the large magnetic Reynolds 
numbers of astrophysics. Nonetheless, large-scale astrophysical 
fields are believed to be predominately generated by so-called 
a-Q. dynamos, in which most of the field amplification occurs 
through the shearing of field lines by ubiq uitous differential ro- 
tation , a process known as the Q. effect (ISteenbeck & Kraus5 . 
Il969l) . For example, many models of the solar dynamo invoke 
the strong shear found in the tachocline at the base of the con- 
vection zone (see, e.g.. ICharbonneaij. l20 1 Q) . Shear alone cannot 
drive dynamo action however, and the a effect, caused by heli- 
cal motions, provides the necessary twist of the sheared field to 
complete the magnetic field amplification cycle. In the Sun, an 
a effect is provided via kinetic helicity due to the interaction of 
stratified convection and solar rotation. 

The a effect can drive a dynamo by itself, being then of the 
so-called type. These dynamos are of great theoretical inter- 
est due to their simplicity, but are expected to be outperformed 
by a-Q dynamos in the wild. Strictly speaking, a-Q. dynamos 
should be referred to as a^-Cl dynamos as the process of 
course continues to occur in reality, even in the presence of the 
Q effect. However, in the mean-field approach one sometimes 
makes the so-called "a - Q " approximation by neglecting the 
production of toroidal field by the a effect entirely in favor of 
the Q. effect. This also applies to the present paper where we 
consider numerical solutions of the original equations in three 
dimensions with turbulent helical flows. However, we will nev- 
ertheless refer to a-Q. and a^- Q regimes when shear is domi- 



nant or comparable with amplification by the helical turbulence, 
respectively. 

Very often, a linear stability analysis of a given setup re- 
veals that several different dynamo modes are expected to be 
excited at the same time. While during the linear stage the rela- 
tive strengths of these modes are determined by the initial con- 
ditions, the mode or mixture of modes of the final saturated state 
is decided by the quintessentially nonlinear interactions between 
the modes in their backreaction on the flow. The naive guess that 
the final state should always be characterized by the mode with 
the highest growth rate , has turned out not to be valid in general. 
In lRadleret dl ([T990l) it was shown for a mean-field dynamo 
model with anisotropic a that within the appropriate parameter 
range both axisymmetric equatorially anti-symmetric and non- 
axisymmetric equatorially symmetric modes can be stable so- 
lutions of the non-linear system. For a system with differential 
rotation it was also shown there that the stable solution can well 
be a mixture of axisymmetric and non-axisymmetric modes. 

In direct numerical simulations of a geodynamo model with 
stress-free boundary conditions, it has been observed that again 
two different dynamo s olutions, a dipo lar and a "hemispherical" 
one, can both be stable dChristensen et al., 1999; Grote & Buss^ 
2000). Because of the free fluid surface in that model, this might 
even be taken as a hint for the possibility of non-unique stable 
sta tes in stellar setups as well. 

iFuchs etaP (119991) have demonstrated an even more extreme 
case with a dynamo powered by a forced laminar flow. In the 
course of the magnetic field growth, the Lorentz force arranges 
the flow into a different pattern, which is hydrodynamically sta- 
ble, but unable to drive a dynamo. As the dynamo dies out sub- 
sequently without a chance to recover, it was named "suicidal". 

Hence, the question for the character of the final, saturated 
stage of a dynamo cannot reliably be answered on the basis of 
a linear approach and the study of the nonlinear model might 
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unveil very unexpected results. Here, we will show in a sim- 
ple setup that, while a-Q dynamos do grow faster than dy- 
namos, non-linear effects are capable of driving transitions from 
a -f2 modes to modes. As the two competing dynamo modes 
are excited for the same parameter set, i.e., the same eigenvalue 
problem, we refer to them as fratricidal, in reminiscence of the 
aforementioned suicidal dynamos. 

The two astrophysical dynamos for which we have long 
time-series, the solar dynamo and that of the Earth, both ex- 
hibit large fluctuations. The solar dynamo in particular is known 
to go through pr olonge d quiescent phases such as the Maunder 
minimum (Eddv, 1976). A conceivable connection with fratrici- 
dal dynamos makes understanding how non-linear effects define 
large-scale dynamo magnetic field strengths and geometries a 
matter of more than intellectual curiosity. 

In Section |2] we sketch the mean-field theory of and a^- 
Q dynamos. In Section [3] we describe our numerical set-up and 
briefly discuss the test-field method, a technique to extract the 
turbulent transport coefficients of mean-field theory from direct 
numerical simulations. In Sections|4]and|5]we describe different 
transition types, and we conclude in Section |6] 

2. Mean field modeling 

In the magneto-hydrodynamic approximation, the evolution of 
magnetic fields is controlled by the induction equation 



^ = Vx(t/xB-77/), 

ot 



(1) 



where B is the magnetic field, J - V x B is the current den- 
sity in units where the vacuum permeability is unity, and 77 is the 
microphysical resistivity. A common approach to dl} is mean- 
field theory, under which physical quantities (upper case) are 
decomposed into mean (overbars) and fluctuating (lower case) 
constituents: _ 

B^B + b. (2) 

The mean can be any which obeys the Reynolds averaging rules, 
and is frequently assumed to be a spatial one filtering out large 
length-scales (a two-scale approach). Here we will however use 
planar averaging, either over the xy plane so that B - (B),iy - 

B(z) = B^ or over the yz plane, that is, B - (B),., = B(x) = B^, 
where (■)^ denotes averaging over all values of the variable ^ (or 
volume, if not specified). A mean defined by averaging over y 
only will also be used. 

Under Reynolds averaging Eq.O becomes 

^^VxiUxB + E-Tjll (3) 
ot 

— = Vx(Uxb+uxB + uxb-6- T]j), (4) 
ot 



where £ = « x ft is the mean electromotive force (EMF) associ- 
ated with corTelations of the fluctuating fields. 

Symmetry considerations allow one to write the £ as a func- 
tion of the mean-fields in the system. In the case of a planar 
averaging scheme, the equation becomes 



&i = aijBj-j]ijJj + ' 



(5) 



where aij and 77,^ are turbulent transport coefficients, and aver- 
aged quantities depend on one spatial coordinate only. The tra- 
ditional a effect is described by the symmetric part of the ten- 
sor aij, and requires helicity in the flow. The symmetric part of 



describes turbulent dissipation, and, in the isotropic case, ap- 
pears equivalently to the microphysical resistivity 77. It is there- 
fore termed the turbulent resistivity, rjt- When assuming that £ 
can be completely represented by the mean magnetic field and its 
first spatial derivatives, the Taylor series in (|5) can be truncated 
after the term rjijjj. A more complete formula would include 
higher spatial as well as temporal derivatives. 

2. 1 . Mean-field dynamo action 

Let us assume a large-scale shearing flow of the simple form 

Us = Sxy. (6) 

and velocity fluctuations which are isotropic, homogeneous, and 
statistically stationary. Consequently, if and 77,^ are assumed 

to be independent of B (the kinematic limit), then they reduce to 
constant scalars a and rjt^ 

If this system were to contain a y-dependent mean field, the 
shear would induce field constituents which are proportional to 
X. We restrict ourselves here to periodic spatial dependencies and 
hence exclude such unbounded fields. The evolution of harmonic 
mean magnetic fields is given by the solution of the eigenvalue 
problem 

-rjjk^ -iak- 
iak, + S -rjjk^ -iak^ 
iak I- —rjjk^ 



AB 



B, 



(7) 



where B - B exp(i^ ■ x + At), rij = % + 77, and k^ — k\ + k^. The 
resulting dispersion relation reads 



(A + T]T:k^)[(A + Tjjk^f - a^k^ + iaSk-] 



0, 



(8) 



with eigenvalues (apart from the always decaying modes with 
B, = 0) 

A^ = -T]Tk^ + {a^k^ - iaSk,^^'^ . (9) 

It can easily be seen that there are two "pure" modes with par- 
ticularly simple geometries: the mode with k^ - does not 
depend on S and has the form 



B = B"" (0, sinfe^x, + cos k,x) , 



(10) 



where the growth rate is A"" - \akx\ - rjjkl and B"" is an am- 
plitude factor. The upper (lower) sign corresponds to positive 
(negative) ak^. 

In contrast, the a^- Q. mode with k^ = does depend on S 
and has, for S » ak- (the a-Q. approximation) the form 



B " = B"" {sm[k,(z - ct)], V2 
c - sign(Q'5) -^JlaSjlic^] . 



sin[^j;(z - c, 



f) + 0],o), (11) 
(12) 



In the above B"" is, again, an amplitude factor, represents, 
for S > {S < 0), the +7r/4 (±37r/4) phase shift between 
the X and y components of the mean field, and upper (lower) 
signs apply for positive (negative) v alues of ak,; see Table 3 of 
[Brandenburg & Subramanianl(l2005b . The corresponding growth 

rate is 

^R{r") = ^J\aSk,\/2 - 77T/t? . (13) 



' Strictly speaking, shear could introduce anisotropy felt by mean 
fields with non-vanishing z-components. Our results do not reveal any 
such. 
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Fig. 1. Time series for a dominantly a -D. dynamo with Re^ - 
20, Ptm - 5 and k{ ^ 3.1. Left: rms value of defined as 

Q. mode (black/soUd), 

,2 



{B^ )-' , to be associated with the a^- 

and of B^, defined as {B^ )]!^, to be associated with the a 

— z 

mode (red/dashed). Right: butterfly diagram of B^. showing the 
dynamo wave of the a^-Q. mode. 



For equal |^^|, the a-D. mode grows faster than the mode. 

A key characteristic of a^-Q. solutions is that the growth 
rate A has a non-vanishing imaginary part k,c which results 
in traveling waves with phase speed c. The wave nature of 
a^-Q solutions is a significant draw in explaining the oscilla- 
tory solar magnetic cycle. For a characteristic a^-Q. dynamo 
found numerically with a setup described below, we show in 
Fig. [T]the time-series of rms values of B alongside the travel- 
ing wave in the z-t plane ("butterfly diagram"). This sol ution is 
similar to those considered recently by Kapyla & Bra ndenburg! 
(l2Q09h . There are other sources for such oscillations how- 
ever. Admittance of a spatially varying a enables oscilla- 
tor y and hence traveling wave solutions in pure dynamos , 
seelBarvshnikova & ShukurovYigS?'). Radler & Brauerl(ll987h . 
[Stefani & Gerbeth (2003) , Mitra et al. (2010). 

The mean fields of modes are force free, while a^-Q. 
modes cause a potential force which has minimal effect as long 
as the peak Alfven speed is sub-sonic. Within kinematics, the in- 
duction equation allows for superimposed and a^- Q modes 
and in Section |4] we will discuss the interactions within such a 
superposition. 



3. Model and Methods 

3. 1 . Numerical setup 

We have run simulations of helically forced sheared turbulence 
in homogeneous isothermal triply (shear) periodic cubic do- 
mains with sides of length 27r. The box wavenumber, which is 
also the wavenumber of the observed mean fields, is therefore 
ki = I. Unless otherwise specified, our simulations have 64^ 
grid points. For the shear flow we have taken the one defined by 



^ When assuming both and k. to be different from zero, but keep- 
ing the a -CI approximation valid and k- fixed, the phase speed of the 
dynamo wave does not change while the growth rate is reduced by rjk^^. 
However, the eigenmode has now a z component ~ -k^jk-By. Such 
modes were not observed in our simulations. 



(|6]l. We solve the non-dimensionalized system 
dA 

^UxB + nV^A 

dt ' 

DU , 1 
-— = -c,Vlnp+ - JxB + Fy 
Dt p 



Dlnp 
Dt 



= -V 17, 



(14) 



(15) 



(16) 



where Cj = 1 is the isothermal sound speed, p the density, f vise = 
p"' V ■ (2pvS) the viscous force, = jiUij + Ujj) - ■ U 

is the rate of strain tensor, v is the kinematic viscosity and / the 
forcing term. We use the Pence. Cod^ which employs sixth- 
order explicit finite differences in space and a third order ac- 
curate time stepping method. While our code allows full com- 
pressibility the simulations are onl y weakly com pressible (small 
Mach number). As in earlier work (Brandenburg 2001 ), in each 
time step the forcing function is a snapshot of a circularly polar- 
ized plane wave. All these waves have the same handedness, but 
their direction and phase change randomly from one time step to 
the next. This forcing provides kinetic helicity. The wavevec- 
tors are taken from the set of vectors that satisfy periodicity 
and whose moduli are adequately close to the target forcing 
wavenumber kf. 

The magnetic vector potential is initialized with a weak 
Gaussian random field, the initial velocity is given hy U - Us 
and the initial density is uniform. In Table [1] we have collected 
the control parameters and some key derived quantities of the 
model. Two parameters of note are the magnetic Reynolds and 
Prandtl numbers. 



Re 



M 



^jTjkf, Ptm = v/t]. 



(17) 



To characterize the turbulence, we provide values of a and 7/t 
which characterize the corresponding tensors as described in 
Section |2] These were determined using the test-field method 
with test-field wavevector k = x or k - z- 

For our purposes, we require the helical turbulence to be 
strong enough that the dynamo can safely be excited. For 
this we guaranteed that in all our simulations, Re^ is above 
the critical value (of the order of unity) for dynamos in the 
corresponding shearless setup. Further, some of the transitions 
we will study require long simulation times due to their rar- 
ity, which constrains us to modest numerical resolutions. This 
in turn prevents our (explicit) numerical resistivity from being 
small, so the turbulent velocities must be reasonably large for 
the stated super-critical values of Re^. Choosing furthermore 
subsonic shear speeds, we are restricted to a modest region of 
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Table 1. Control and derived parameters 



V 


Control par. 


Microphysical viscosity 


1 


Control par. 


Microphysical resistivity 


S 


Control par. 


Shear (Us = S xy) 




Control par. 


Forcing amplitude 


k, 


Control par. 


Forcing wavenumber (generally kf « 3.1) 


PrM 


v/ri 


magnetic Prandtl number 






RMS turbulent velocity 




UrmJrjkf 


Magnetic Reynolds number 


k, 


k, = 1 


Wavenumber of mean fields (box wavenumber) 


«res 


l/rik] 


Resistive time (mean fields) 


^turb 


1 /^rms^f 


Turbulent time 
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parameter space. In light of these limitations we operate mostly 
in a PrM > 1 regime. 

3.2. Test-field method 

A fundamental difficulty in extracting the tensors and rjij 
from a numerical simulation of (fT4l i-(fT6Tl is that (|5]l is under- 
determined. Turbulent transport depends on the velocity field, 
so "daughter" simulations of the induction equation, whose ve- 
locity fields are continuously copied from the main run, share 
the same tensors aij and rjij. It is therefore possible to lift the 
degeneracy by running an adequate number of daughter sim- 
ulations with suitably chosen "test" mean fields. We employ 
this test-field method (TFM ); f or an in depth overview see 
ISchrinner et aH (l2005l l2007l) and iBrandenburg et al.l (l2008allbh . 
Recently the original method has been extended to systems 
with rapidly evolving mean-fi elds, requiring a more compli- 
cated ansatz than Eq. alphaeta dHubbard & BrandenburgI l2009l) 
and to the situation with magn etic background turbulence 
dRheinhardt & Brandenburgll20I0l) . 

In addition to calculating planar-averaged turbulent tensors 
as described in the references above, we will be interested in 
tensors that depend both on x and z (that is, are y-averages). For 
this, we generalize (|5j to 

dBj 

6i^aijBj+/3ijk^ + -- - . (18) 

There are 27 tensor components (as dyB - 0), so nine test-fields 
are required, which we choose to be of the form 

Bi'i ^ B'^fq(x,z)6ipei, p= 1,2,3, ^ e {cc, sc, cs), (19) 

where fq(x, z) is defined, according to the choice of q, to be one 
of the following functions: 

coskixcoskiz, sin/cixcos^iz, cos^ixsin^iz, 

and is, as standard for test-field methods, an arbitrary ampli- 
tude factor Although the wavenumber of the test fields is usu- 
ally treated as a varying parameter we need here to consider 
only the single value ki because the fastest growing and also 
the saturated dynamos in the simulations are dominated by this 
wavenumber, the smallest possible in our periodic setup. As is 
often the case in applications of the test-field method, we will 
occasionally be faced with unstable solutions of the test prob- 
lems. We treat that difficulty by pe riodically resetting the test 
solutions (see H ubbard et aUl2009 l). Since it takes a finite time 
for the test solutions to reach their stationary values, and as this 
time is frequently close to the required reset time, only limited 
windows in the time series of the data are valid. 

4. Dynamical interactions of and a^-Q. modes 

Here we report on the results of our simulations a first set of 
which is characterized in Table |2] In Figure |2] we show time se- 
ries for Run A, which saw a transition from a z-varying a^- O 
dynamo (B^) to an x-varying dynamo (B^). As is made clear 
in the bottom panel, there was a prolonged period where the 
two modes were coexisting while their relative strengths were 
changing monotonically. However, note that B^ is stronger than 
B^, that is, the field is distorted during the transition. Run 
A was repeated 16 times with the same parameters, but differ- 
ent random seeds, and all these runs exhibited similar behavior. 



fratricide of a ~ CI dynamos by their siblings 

Mean magnetic fields B^, y — component 




Fig. 2. Time series for Run A. Upper row: same quantities as 
in Fig. [1] Middle row: B^. and Z?, , to be associated with the 
mode. Note that the a^- Q and modes coexist during the 
transition. Lower panel: rms values of the components of B^ and 
B^. 



Likewise we performed runs where both the value of and the 
numerical resolution (cf. Runs B-D, I,J) were varied. As these 
additional runs also showed the same transition pattern, we con- 
clude that it is deterministic for this level of shear and forcing. 
More, we conclude that for these cases the a^- Q mode is unsta- 
ble to the growth of an mode due to non-linear effects. Runs 
with the dynamical parameters (S , Mrms) of Table |2] inevitably 
generate fields from a^-Q. fields after modest times, runs 
with significantly different parameters will usually (for most of 
the random seeds) exit the kinematic regime into an a^- Q mode, 
and stay in that mode for a prolonged time with no sign of an 
field. Nonetheless even such simulations can occasionally fail 
to fully enter in the a^- Q regime, instead exiting the kinematic 
regime into an mode, as shown in Figure [3] 
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0.0 0.5 1.0 1.5 2.0 
t r) k? 



B , y-component 



0.5 1.0 1.5 2.0 
t r, k? 



-0.3 0.0 



B**, z-component 





0.5 1.0 1.5 2.0 
t 7J k? 



-0.2 0.0 0.2 



0.5 1.0 1.5 2.0 
t ri k? 



-0.2 0.0 0.2 



Fig. 3. Time series for Run E that never entered a quasi- 
stationary a^-Q. regime. Top left: rms values of to be asso- 
ciated with the Q-^-Q mode (black) and of to be associated 
with the mode (red). Note the considerably faster growth of 

the a^- Q mode during the kinematic phase. Top right: Butterfly 
— z 

diagram of B^, , showing the traveling dynamo wave during the 
kinematic phase only, but later merely fluctuations. Bottom: 
components of B^. 



4.1. Mean-field approach 

Clearly, the transition from an a^- Q. mode to an one must be 
a consequence of the back-reaction of B onto the flow. Within 
the mean-field picture, there are two channels available for it: 

(i) the back-reaction onto the fluctuating flow, usually described 
as a dependence of ajj (more seldom rjjj) on the mean field and 

(ii) the back-reaction onto to the mean flow by the mean Lorentz 
force, which might again be decomposed into a part resulting 
from the fluctuating field, j x b, and one resulting from the mean 



Table 2. Run parameters 



Run 


Res. 


-5 




-a' 


nl 


Rsm 


PrM 


T* 


Run A 


64' 


0.05 


0.11 


0.04 


0.025 


37 


5 


2-3 


RunB 


643 


0.05 


0.17 


0.04 


0.03 


26 


2.5 


2 


RunC 




0.05 


0.14 


0.04 


0.027 


44 


3 


4 


RunD 




0.05 


0.14 


0.04 


0.027 


90 


6 


1.5 


RunE 


643 


0.02 


0.13 


0.04 


0.023 


90 


10 


N/A 


Run I 


643 


0.05 


0.15 


0.04 


0.036 


49 


1 


3 


Run J 


643 


0.05 


0.19 


0.04 


0.035 


31 


0.5 


1.5 



Notes, t Time-averaged values determined through the TFM using har- 
monic test fields with k = x oi k = z- The results are identical due to 
homogeneity of the time-averaged turbulent velocity, i r = fdur/fics is 
the duration of the transition of the type discussed in Section|4j counting 
from multiple random seeds for Run A. 



field, J X B. Here, we will deal with a flow generated by the 
latter force that straddles the distinction of means and fluctua- 
tions: it survives under y-averaging, but vanishes under the xy 
and yz averaging that reveals the a^- Q. and dynamos respec- 
tively. For simplicity we consider magnetic field configurations 
that would result from a superposition of linear modes of the 
ff-Q and dynamos, given in equations (fTTI) and ( fTOb respec- 
tively. Such a situation will inevitably occur during the kinematic 
growth phase if both dynamos are supercritical, but is only rele- 
vant for analyzing the back-reaction onto the flow if it at least to 
some extent continues into the non-linear regime. Our analysis 
is linear in nature, so while it provides a qualitative framework 
for understanding the transition process, it is surely not quanti- 
tatively accurate. 

In order to be able to consider both B ^ and B ^ as mean fields 
under one and the same averaging, we have now to resort to y 
averaging. Moreover, for the sake of clarity we will occasionally 
subject the resulting x and z dependent mean fields further to 
spectral filtering with respect to these coordinates. That is, we 
will consider only their first harmonics ~ e'*'*'^~' as mean fields. 

Let us represent the mean field {B)y as superposition of a 
B^ resembling the (x varying) mode B"" (Eq. (fTOl )) and a B^ 
resembling the (z varying) a-Q mode B"^ (Eq. (fTTTl): 



B^ = 



B^ = B^ 



sinikiz') ^ 
G sin(feiz' + (j)) 




Hsmk\x 

, cos kyX , 



-G cos(feiz' + (p)\ 
cos(^iz') 




^ 

sinA:ix 
,//cos kix, 



(20) 



with zf = z - ct recalling that c is the speed of the dynamo wave 
(Eq. ((nil). In the above, 7t/4 < (p < 37r/4 and G, H,ki > 
are appropriate for a > 0, The parameters G and H capture the 
difference in the strengths of the y and z components (or^ ) or the 
X or y components {a^- Q ), respectively. We expect G > 1 as 
shear amplifies the y component of an a^- Q mode well above its 
X component. The inclusion of the parameter H, which is unity 
for pure modes will be justified below, but can already be seen 
in the different strengths shown in Figure |2] lower panel. 

The mean Lorentz force J x B for the superimposed fields 
can be written as 

Fl = {FL}y = (21) 

kiB^B^ cos kix[G cos(kiz + 4>) + H ^m{kiz!)]y + V<D. 

As the Mach numbers were found to be small throughout, we as- 
sume incompressibility and hence drop the potential component 
VO. Further, we assume that Fi and the mean velocity driven 
by it are simply linked by a coefficient K ^ l/vjk^, where the 
total viscosity vt is the sum of the molecular v, and the turbulent 
viscosity Vt. Thus we can approximate the mean velocity due to 
the interaction of the superimposed mean fields as 

Ul - Ul cos kix [G cos(kiz' + (p) + H sin(A:iz')] (22) 

where Ul - Kk\B^B^. Clearly, this flow, having merely a y 
component, shows quadrupolar geometry in the x - z plane as 
Ui.y can be rewritten in the form t/^ cos kixcosikiz' + 4>') with 
a new amplitude and phase, U'^ and (/>'. 

The simulations show indeed a dominant part of that shape 
in the Lorentz-force generated mean flow as can be seen from 
Figure |4] There the y-averaged Uy is shown together with its 
Fourier constituent ~ e'*^'^*"^""'. The latter contains approximately 
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y-Qveroged U 

We can write 

Hcos(kiz) - G sin(feiz' + </.)= (27) 
{{H cos <p - G) sin(A;iz' + <^)] i + {H cos ^ cos{k\z' +</>)] 2 ■ 

If //cos (p - G < 0, as expected since // ~ 1, G > 1, term []i 
in (IZTT i will act to damp Bf, that is, the perturbative a^- Q wave. 
Further, term []2 is opposite in sign to the time-derivative of such 
a wave, so it slows or reverses the direction of wave-propagation. 

4.4. Mean-Field Evolution 

Here we assume again domination of the a^- Q. mode, that is, 
» B^. With Eqs. ^ and (l24l) the eigenvalue problem for the 
modified field is then (adopting - ki, k. - 0) 




Fig. 4. {Uy)y for Run A, taken at eai'ly time (t = 1.45fres) when 
is still of only modest strength. Plotting area is shifted in 
X and z to make the quadrupolar geometry clear. Overplotted 
contours; quadrupolar constituent ~ coskix' coskiz' 

10% of the energy in this component, or U 'j - J/y,rms/3, indi- 
cating that the assumptions made in deriving (l22b are reasonably 
well justified in a non-linear system. 

Upon interaction with a B'^ or a of the form (|20] |. the 
mean flow Ul in (l22l l generates an S^iz) and £j(x), respectively. 

4.2. Dominating a^-£l mode 

If B^ » B^, then can be treated as a perturbation, and we 
can drop higher order terms in B^. Accordingly, the z-averaged 
EMF due to the flow Uiis 



£^ = (UlxB^), = ^^B^B^^iGsiiKp- H)coskix t 



(23) 



The curlof fliisEMFis 



V x£^ = B'^Ism(kix)y, 



Kk^B^^ 
I = — \- (Gsin^ 



•//). (24) 



If G sin > H, then / > and for // > this EMF reinforces 
By - B^H sin kix. Thus we see that the inclusion of the param- 
eter H in the ansatz for B^, Eq. ( l20l l. was justified as By receives 
enhanced forcing in comparison to B?. 

4.3. Dominating mode 

If B^ » B^ then we can in turn treat B^ as a perturbation. 
Further, as the system is dominated by the mode, we will 
have H ~ 1 . In this case we find 



(25) 



= -^B^B^ [G cos(kiz' + (p) + H sin(;fciz')] x. 



and 



Kk^ 

V X £^ = —^B^^ [Hcosikiz) - G sm(kiz' + 0)] J. (26) 



A^B^ 



-rirk] 
S 





iak\ 





-i(aki + 1) 
-riTk\ 



B' 



with eigenvalues 



= -rik\ + ^Jakl(akl + I). 



(28) 



(29) 



Making the approximation / » ak[, similar to the a-Q. approxi- 
mation S » aki , we find 



-rjkl + -\JaIki 



(30) 



The above should be compared with the growth rate of the a-Q. 
dynamo, A"^ from (fT3l l which is not touched by the occurrence 
of /. The a-Q. dynamo saturates when a has been quenched 
such that the product aS settles at the marginal value la^ | = 
2Tij\ki\-^. If the parameter / becomes comparable with the shear, 

i.e., I ~ S , then B^ might grow even when the a^-Q field is 
saturated, i.e. A^ > %(A"^) - 0. In other terms, the a--Q mode 
is unstable to the growth of ?l fratricidal a^ field, so the transition 
will take a well defined time from the onset of the non-linear 
stage, determined by A^ . 

We test this theory for Run A at the time of Fig. |4] f = 
1.45f,es, extracting G and H from the relative strengths of the 
X and y or y and z components of the averaged fields B^ or 
B^, respectively, after a projection onto the first harmonics; see 
Eq. (I2OI 1. The parameter / is calculated from the magnetic and 
velocity fields using 



2B^ 



(G sin0-//), 



(31) 



wifli Ul = Uy ^JG^ + -2GH sin (p, where t/[ is the ampli- 
tude of the quadrupolar constituent of the velocity field seen in 
Fig. H We find U'^ ^ 0.07, H ^ 2.9, G ^ 4.9, / ^ 0.09, and 
confirm that ^ n/A. As I > S - 0.05, the growth of the 
jc- varying mode even when the a-Q mode is saturated is not 
surprising. Repeating this run (keeping the control parameters 
fixed) 16 times with different random seeds changed the occur- 
rence time of the transition by only one resistive time, suggesting 
that the transition is an essentially deterministic process. 

We have never seen a reverse transition from the a^ state 
back to the a^-Q state. This may be understood in terms of in- 
teracting modes, with the a^- Q mode being suppressed once the 
a^ mode is dominating; see Sec. 14.31 
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1 L 



4 6 8 10 12 

t r, k? 

Fig. 5. Time series for Run F (solid lines), with rms values of 
(black) and (red). Broken lines: restarts from the main 
run with new random seeds, vertically offset for visibility. All 
the runs end up with the same energies in B^ and B^. Vertical 
lines: restart times. 

Mean magnetic fields 




t r) kt 

Fig. 6. Time series of Run G (solid) along with a sibling run 
(dashed) with different seeds, showing significant differences in 
the transition start time. Dash-dotted: a run which never entered 
the a^- regime. 



Table 3. Run parameters 



Run 


Res. 


-S 




-a+ 


nt 


Rsm 






RunF 




0.2 


0.14 


0.02 


0.1 


90 


10 


5,9* 


RunG 


643 


0.2 


0.099 


0.01 


0.16 


63 


10 


5, 20* 


RunH 


643 


0.1 


0.085 


0.017 


0.037 


27 


5 


25 



Notes, t see Tabled $ t = ftrms/'res is the time when transition occurred; 
for F and G over multiple realizations with differing random seeds. See 

Fig.m 



5. Random transitions 

Not all transitions fit the above deterministic picture of inter- 
acting a^-Q and modes. In Fig|5]we present a set of time 
series of the rms values of B^ and B^, all related to Run F of 
Table [3] Secondary runs were performed by branching off from 
the original simulation either at f = 5fres, when the a^-Q. mode 
is well established and stationary, or at f = 8.25fres, immediately 
before the transition to the mode is launched. The only dif- 
ference between all these runs is in the random seed, which is 
used by the forcing algorithm. In all, the time until the transi- 
tion starts varies by x 2.5f,es, and many more turbulent turnover 
times (Re^^f /^i - 800 turbulent turnover times per resistive 
time). The time elapsed during a transition is always of the or- 
der of fies/2, unlike 3fies for the process seen in Fig. |2] Thus it 
is suggestive to assume that there might be a very slow, still es- 
sentially deterministic process, preparing the transition, which 
is likely resistive in nature as that is the longest obvious "na- 
tive" timescale of the system. Slow resistive effects are known 
to exist in dynamos, for example the slow resistive growth of 
dynamos in periodic systems. However, transitions can in- 
deed occur at very different times including the extreme case in 
which a run never develops a quasi-stationary a^- Q mode, but 
instead enters the state almost immediately after the end of 
the kinematic phase, see Fig. |6] (run G of Table O. We believe 
therefore that under certain circumstances the transition process 
is not a deterministic one, in that it is impossible to predict or at 
least estimate the time until the transition. FigureQis a synopsis 
of simulations that belong to that type, hence do not show the 
instability discussed in Section H) Note that, while correspond- 
ing setups without shear are known to enable modes for the 
entire parameter range, the mode is possibly sub-critical for 
Rcm = 10,5 = -0.1. 

This is different from the interacting mode picture of Sec. 
in several interesting ways. Firstly, the a^-Q mode is here at 
least meta-stable against growth of the mode, as evinced by its 
prolonged life-time (hundreds of turbulent times) and the small 
magnitude of B^, which further is not dominated by a mode. 
A reasonable working hypothesis for the cases of Sec. |4]is then 
that, there, the mode is the only stable solution and, as soon as 
the nonlinear stage has been entered, it starts to devour the a-- 
Q one, settling after a time which is related to basic parameters 
of the system and hence not random. In contrast, for the cases 
considered here, we conclude that both the and the a^- Q so- 
lutions are indeed stable (not only metastable) and the latter has 
a well extended basin of entrainment. Due to its higher growth 
rate the system settles first in the a^- O mode and suppresses the 

mode efficiently. A transition to the latter can only occur if 
a random fluctuation in the forcing is strong enough to push the 
system over the separatrix into the basin of entrainment of the 
or^ mode. This can happen after a rather long time only or im- 
mediately after the end of the linear stage which has both been 
observed. 

Given that the examples for the first scenario (Table |2) dif- 
fer from those for the second (Table |3) mainly in their lower 
rate of shear, our conclusion seems reasonable as stronger shear 
should result in a clearer preference of the a-- D. mode as the 

mode does not feel the shear. Or, in other terms, from a cer- 
tain shear rate S on, the a^- Q mode should acquire a basin of 
entrainment with a finite "volume" that grows with S . If this pic- 
ture is true, transitions in the two scenarios should have clearly 
different characteristics, and indeed, the transition in Fig. |5]is 
markedly faster than that seen in Fig. |2] 
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Fig. 7. Synopsis of runs which did not exhibit the instability 
discussed in Section |4] Runs at the same position differ only in 
random seeds. Circles/black: a significant a^- Q. mode was never 
developed (cf. Fig. [3), size indicates the corresponding number 
of runs (1, 2 or 3). Square/blue: a transition occurred, size rep- 
resents the time until transition (4 to 25f,es)- Diamond/red: the 
a^- O stage was entered, but no transition occuiTed. Size repre- 
sents the time span of simulation (5 to 35fi-es). 



As in the transitions discussed in Sec. |4] we have not here 
seen the mode transit back into the a^-O mode. Some at- 
tempts were made to provoke this reverse transition by perturb- 
ing the state with a (sufficiently strong) a^- Q. mode. While in 
some runs it indeed took over, velocities were attained for which 
the numerics are unreliable, and often proved numerically unsta- 
ble, making the results inconclusive. However, such a behavior 
is not entirely surprising as the a^- Q. saturation process can any- 
way be somewhat wild, cf. Fig. [3] 

The absence of spontaneous reverse transitions appears plau- 
sible insofar the time variability of the mode is much smaller 
than that of the a--D. mode, which can clearly be seen in Fig. [8] 
for Run H. That is, events capable of pushing the system over 
the separatrix are simply much rarer. Significantly longer inte- 
gration times are likely needed for their eventual detection, but 
it is also conceivable that the triggering event never shows up. 

5.1. Large scale patterns 

Run H will be examined here in more detail. Curiously, (f/,), 
taken just during the transition as shown in Fig.|9]does not show 
the quadiTipolar pattern of Fig. 2] It is therefore not surprising 
that the butterfly diagrams in Fig. [10] do not show a direct tran- 
sition from the a^- Q. to the dynamo, as B? develops signifi- 
cantly later than By. This is clearly visible in Fig. |8] lower panel. 
As consideration of the mean flow due to the Lorentz force of the 
mean field alone is obviously not fruitful in explaining this tran- 
sition, we recall that the back-reaction of the mean field onto the 
turbulence opens another channel of nonlinear interaction. 

According to elementary mean-field dynamo theory, the a 
effect is caused by the helicity in the flow: a ~ {w ■ u), where 
w s V X H is the fluctuating vorticity. Further, the back-reaction 
of the mean field on the turbulence, which saturates the dynamo. 




Field components 




25.0 25.5 26.0 

t rj 

Fig. 8. Time series of Run H. Upper panel: rms values of and 
B^. Note the long time before the transition starts in comparison 
to Run F (see Fig.|5]l and the dramatic difference in the fluctua- 
tion levels before and after the transition. Lower panel: rms val- 

— z 

ues of the components. Note the strong difference between By 

— z , 
and Bj. , expected for an a -Q field. More significantly, notice 
^ 

that By develops before B_ . 



is assumed to be captured by the cuiTent helicity (j ■ b). It is 
often related to the magnetic helicity {a b) and thought to reduce 
the original a by producing a magnetic contribution of opposite 
sign. In Fig.[TT]we present time-series of the power spectra of 
these helicity correlators across the transition. We see no clear 
signal around the tx^ansition event. 

5.2. Mean-field modeling with y averaging 

To examine the problem more closely, we recall Eq. [T8]for when 
the mean is defined by a y average: 

£>i(x, z) = aij{x, z)B j(x, z) + jSijkix, z)B j^kix, z). (32) 
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t = 25.2 



M 




Fig. 9. {Uy}y taken during the transition of Run H shown in 
Fig. [8](f = 25.2fres)- Note the lack of a quadrupolar geometry. 
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B*', Y-component 



25.0 25.5 26.0 



B , Z-component 



X 




25.0 25.5 26.0 
t 7j 



25.0 25.5 26.0 
t 7j k? 



Fig. 10. Butterfly diagrams for Run H, (see Fig.|8). Note that B,, 

—X — , 
develops before , i.e., B does not transit from an a -Q field 

straight to an one. 



It is clear that the Fourier constituents of ajj and /Sijk with 
wavenumber in both x and z (the quadrupolar constituents) 
can create an emf out of a field B^, both with the same 
wavenumber: 



(33) 



where the superscripts indicate the coefficients to be the Fourier 
constituents ~ e'*^'**^^' and B? is assumed to vanish. Note that 
each of them is actually given by four values, e.g., the two am- 
plitudes and phases in: 

a-j = a^jj cos(feix + 0;^) cos kiz + a]j cos{k\x + <p]j) sinfeiz. (34) 



Power (A • B) 



Power (J • B) 



Power (W 




Fig. 11. Time series of the helicity power spectra for Run H. 
Horizontal line: forcing wavenumber kf ^ 3.1. Vertical line: 
border between low-resolution observations (every O.Sfies) for 
t < 25fies and higher-resolution observations (every O.lfres ) for 
t > 25fres- Possible features in the latter range are likely due to 
the increased temporal resolution. 



The coefficients relevant for the generation of By (from £f 
only) are 0-3 [, a\ l,, jg ']-, and /5]^23- We have used the test-field 
method (see Sec. 13.2b to find them and present the results in 
Fig. [12] They turn out to be surprisingly large, when compared 
to the rms velocity (e.g., (a'^^^ + a-^^y^^ > 4Mi-ins) and some may 
show a trend across the transition from the a-Q.to the mode. 
This overall trend is hypothesized to be due to the increase in 
Mrms that accompanies the transition from a stronger a - Q field 
to a weaker field with less potential to inhibit the flow. It is 
interesting that with the exception of a^^, the large transport co- 
efficients are all those which generate an £ out of the f-directed 
field, i.e. out of a field that feels the effect of shear. We speculate 
that these coefficients, with, themselves, explicit A:-dependence, 
feel the shear quite strongly. 



6. Discussion and conclusions 

We have demonstrated that, while a^-Q modes are kinemati- 
cally preferred to modes in homogeneous systems that sup- 
port both, the mode acts in a fratricidal manner against the 
former after the nonlinear stage has been reached. This transi- 
tion can occur in at least two different fashions. Further, we have 
not observed the reverse process. One of the two transition pro- 
cesses, based on superposed a^-Q. and modes, operates in 
a basically deterministic fashion through a large-scale velocity 
pattern generated by the interaction of the modes. In contrast, 
we interpret the mechanism of the second process, which may 
start only many resistive times past the saturation of the a^- Q. 
dynamo, by assuming that both the a^- Q. and the modes are 
stable solutions of the nonlinear system. Transitions occur if due 
to the random forcing a sufficiently strong perturbation builds 
up which tosses the system out of the basin of entrainment of 
the a^- Q mode into that of the mode. This hypothesis is bol- 
stered by both the random timing of these transitions and by the 
large time-variability seen in the amplitude of the a^- Q. field. A 
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25.0 25.5 26.0 
t 7) k? 



25.0 25.5 26.0 
t Tj 



Fig. 12. Run H. Upper left panel: rms values of B^, and 
« (Mims - {u^y^^), cf. Fig. [8] Remaining panels: Selected 
quadrupolar moments of aij and j3ijk determined by the test-field 
method and given by ia'if + a]f)^l^ (see Eq. dH)), likewise for 
jSijk- Normalization is by the temporally averaged Mn^s, as Mn^s 
undergoes a slow, steady drift over time. Vertical lines mark the 
times of events visible in the first panel. 



return seems to be much less likely as the level of fluctuations of 
the mode is, by contrast, greatly reduced. 

These results fit with earlier work studying dynamos whose 
non-l inear nature i s fund amentally different from their linear one 
(e.g., iFuchs et all 11999 ). While our simulations are limited to 
Cartesian, cubic, shearing-periodic domains, they are particu- 
larly exciting given that the only dynamo which has been ob- 
sei"ved over a long baseline and which could be either a - or 
a^, the solar dynamo, indeed shows differing modes of operation 
(regular cycles vs. deep minima). The results are also disturb- 
ing in that we have evidence for non-deterministic, rare (as they 
occur in scales of multiple resistive times or hundreds of turbu- 
lent turnovers) mode changes that show no evidence for a return. 
Given that the mode in our simulations seems much calmer 
than the a^- Q mode, a rare random excursion in the field geom- 
etry is likely to be the initiating agent of the transition. While a 
bifurcation between different stable mode s has long been an ac- 
knowledged po ssibility for dynamos (Bra ndenburg et all Il989t 
lJenningslll99Th . a rare, stochastic, possibly uni-directional tran- 
sition is perhaps the most troublesome form of such bifurcations 
except for the ultimate self-extinction. 

The a-Q. dynamo is believed to be common and impor- 
tant for systems like the Sun or accretion disks, which all have 
long life-times compared to turbulent turnover times. It is then a 
daunting possibility that we could be forced to stretch our simu- 
lations over very long temporal base-lines to find the actual long- 
lasting field configuration. More positively, our result, while in 
a different geometry, increases the importance of recent work 
on non-oscillatory a-Q an d oscillatory modes in spherical 
shells for the solar dynamo (iMitra et al.l 1201 Ol ISchrinner et al.L 
I20T1I) . 
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